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Numerical  and  analytic  solutions  in  the  ray-optics  region  have  been  obtained  for  the 
scattering  of  HF  radio  waves  by  a  spherical  electron  cloud  with  a  Gaussian  distribution  of 
electron  density  in  the  radial  direction.  The  results  show  that  a  highly  overdense  or  hard  cloud 
has  a  wide  range  of  backscatter  angles  for  which  the  cross  section  is  almost  constant  with  a 
value  close  to  that  of  an  equivalent  conducting  sphere  (critical  radius  size).  A  slightly  overdense 
or  soft  sphere  has  a  much  narrower  range  of  almost  constant  cross  section  whose  value  is  pro¬ 
portional  to  the  fourth  power  of  the  critical  radius  and  considerably  below  that  of  an  equivalent 
conducting  sphere.  In  the  region  of  forward  scatter,  all  sphere:  have  essentially  the  same  cross 
section,  independent  of  hardness.  For  underdense  spheres,  the  cross  section  is  generally  the 
same  as  that  characteristic  of  forward  scattering  except  in  the  region  of  maximum  deflection 
where  the  cross  section  increases  very  sharply.  It  is  found  that  when  an  underdense  sphere 
becomes  critical,  the  maximum  deflection  is  90°,  whereas  an  overdense  sphere  retains  its 
maximum  deflection  of  180°  at  the  critical  density.  For  overdense  spheres,  the  analytic  results 
are  in  good  agreement  with  those  obtained  numerically;  but  for  underdense  spheres,  they  are 
not  in  good  agreement  at  the  higher  electron  densities. 


1.  INTRODUCTION 

In  recent  years,  studies  of  the  upper  atmosphere 
have  been  conducted  with  artificial  ionized  clouds 
obtained  by  the  release  of  barium  vapor  at  high  alti¬ 
tudes.  To  help  in  the  understanding  and  interpreta¬ 
tion  of  experimental  data,  a  theoretical  study  is  being 
made  of  the  scattering  of  radio  waves  by  an  artificial 
electron  cloud. 

Because  the  electron  cloud,  which  is  several 
kilometers  in  extent,  is  considerably  larger  than  a 
typical  wavelength  in  the  HF  band,  the  method  of 
geometrical  optics  has  been  utilized.  For  those  fre¬ 
quencies  at  which  the  effect  of  the  magnetic  field  is 
not  negligible,  a  detailed  ray-tracing  technique  such 
as  that  furnished  by  the  Haselgrove  method  [Hasel- 
grove,  1955]  must  be  employed.  If,  however,  the 
magnetic  field  can  be  neglected,  the  scattering  may 
be  computed  directly,  at  least  for  the  case  of  spheri¬ 
cal  symmetry,  with  increased  accuracy  and  consid¬ 
erable  reduction  in  computational  labor. 

Our  present  investigation  is  concerned  with  the 
frequencies  at  which  the  magnetic  field  is  negligible. 
We  shall  assume  that  the  distribution  of  electron 
density  is  spherically  symmetric  and,  as  indicated  by 
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diffusion  theory,  has  a  Gaussian  form  in  the  radial 
direction.  The  collision  frequency  within  the  cloud 
is  generally  small  compared  to  the  incident  fre¬ 
quency;  thus  absorption  may  be  neglected  in  our 
calculations. 

2.  ANALYSIS 

A  ray  initially  a  distance  b  from  the  z  axis  is 
deflected  through  an  angle  p  at  a  given  position  r, 
<#>  within  the  cloud,  and  is  scattered  through  an  angle 
0  as  it  emerges  from  the  cloud  (Figure  1).  For  a 
spherically  symmetric  plasma,  the  scattering  is  inde¬ 
pendent  of  meridian  angle  and  can  be  described  in 
terms  of  the  scattering  angle  6  in  the  plane  of  in¬ 
cidence  and  the  impact  parameter  b.  The  relation  be¬ 
tween  b  and  0  may  be  obtained  from  the  scattering 
integral  for  the  two-body  central  force  problem 
[Goldstein,  1 950a]  and  the  analogy  between  particle 
dynamics  and  ray  optics  [Goldstein,  195Ub],  or  by 
direct  consideration  of  the  scattering  of  a  ray  by  a 
spherically  symmetric  plasma  [Kelso,  1964].  We  may 
express  this  relation  in  the  form 

2  b  l,  0*V  -  b Vj175  (,) 
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Fig.  1.  Sketch  showing  ray  path  and  scattering  parametcis. 


where  p  is  the  index  of  refraction,  r  the  radial  posi¬ 
tion  of  the  ray,  and  r,  the  turning  point  or  reflection 
point,  that  is,  the  location  at  which  the  ray  starts  to 
move  away  from  the  center  of  the  sphere.  For  a 
Gaussian  distribution,  the  electron  density  A'  and  the 
index  of  refraction  have  the  form 

X/  Am.,  =  exp  (-//a2)  (2) 

and 


p  -  I  —  A'/ A'.,  -  I  —  exp  [(r,.,2  -  rJ)/aJ\  (3) 

where  a  is  the  Gaussian  radius  and  the  subscripts 
max  and  cr  refer  to  maximum  and  critical  values, 
respectively.  It  is  convenient  to  change  the  variable 
of  integration  to  jc  =  a/r  and  to  take  a  as  the  unit  of 
length  to  yield 


and 


=  b 


r  -  j* 
Jo  (/  - 


(4) 


of  r  so  that  p:  varies  very  rapidly  over  a  small  range 
of  .r  well  inside  the  inflection  point  which  occurs  at 
.t  =  (2/3 ) We  do  not  c.vpcct,  therefore,  the  de¬ 
tails  of  the  variation  to  be  very  important  and  should 
obtain  a  reasonably  accurate  solution.  However,  for 
a  slightly  overdense  plasma  (small  r,.r)  at  small  im¬ 
pact  parameters,  the  turning  point  occurs  at  a  small 
value  of  r  with  a  consequent  slow  variation  of  p?  over 
a  wide  range  of  .r  which  includes  the  inflection  point. 
The  manner  of  variation  is  important  here,  and  we 
anticipate  more  difficulty  and  less  act  uracy  in  obtain¬ 
ing  approximate  solutions.  Since  a  ray  penetrates 
very  little  into  a  highly  overdense  plasma  and  deeply 
into  a  slightly  overdense  plasma,  it  is  convenient  to 
refer  to  these  cases  as  ‘hard’  and  ‘soft’  spheres,  re¬ 
spectively.  \k  hen  the  plasma  becomes  underdense, 
the  turning  point  is,  as  in  the  slightly  overdense  case, 
small  at  small  impact  parameters;  and,  in  addition, 
the  scattering  is  now  confined  to  the  forward  region 
(0  ^  90“ ).  We  may  therefore  consider  the  under¬ 
dense  case  as  a  ‘very  soft’  sphere  and  expect  diffi¬ 
culties  in  obtaining  an  approximate  solution. 

3.  ANALYTIC  PROCEDURE 

In  view  of  the  somewhat  different  treatments  re¬ 
quired,  the  overdense  (hard,  intermediate,  and  soft) 
and  the  underdense  spheres  are  treated  separately. 

3.1.  Hard  sphere.  We  consider  first  a  typical 
hard  sphere  with  r,.T  -  2,  for  which  a  plot  of  /»■ 
against  .r  is  given  in  Figure  2.  Because  of  the  rapid 
drop  in  /r,  it  is  more  accurate  to  use  a  single 
quadratic  term  than  a  combination  of  linear  and 
quadratic  terms  and  write 


=■  I  —  exp  (at,,  2  -  .v"2)  (5) 

The  differential  cross  section  in  the  direction  0,  <r  ( 0 ) , 
is  then  calculated  by 

ff(0)  =  (4ir6/sin  6)(db/d9)  (6) 

Equations  4  and  6  can  always  be  solved  numeri¬ 
cally  for  a  given  variation  of  /  with  x.  An  approxi¬ 
mate  analytic  solution  may  be  obtained,  however,  if 
the  function  p*  may  be  fitted  with  reasonable  ac¬ 
curacy  by  a  linear  or  quadratic  function  of  x,  since 
the  integral  in  (4)  can  then  be  expressed  in  terms  of 
elementary  functions.  The  properties  of  the  cross 
section  and  its  dependence  upon  the  governing 
parameters  in  a  given  region  of  interest  can  then  be 
clearly  exhibited.  For  a  highly  overdense  plasma 
(large  r,r),  the  turning  point  occurs  at  a  large  value 


a  =  I  -  aA.x  -  x,  Y  x  £  Xj 

where  .va  is  the  point  at  which  the  parabola  starts  and 
is  determined  by  noting  where  p*  for  the  Gaussian 
has  not  yet  decreased  significantly  from  unity.  The 
constant  a3  may  be  evaluated  by  noting  that  p  —  0 
at  .vir  =  0.5  to  give  at  =  44.5.  A  plot  of  p1  for  the 
parabola  shows  that  it  is  in  good  agreement  with  the 
exact  curve  over  the  entire  range. 

The  scattering  integral  can  then  be  written  as 


+ 6  f  n^*r--c-sf  m 
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Fig.  2.  Plot  of  ti-  for  hard  sphere  (r„  =  2)  and  approxi 
mation  by  parabola. 


where  A  --  1  -  b-x22,  B  =  2b-x2,  C  -  a  .  +-  and 
u  =  x  -  ,r2  is  the  variable  of  integration. 

Performing  the  integration  and  noting  that  the 
quadratic  vanishes  at  the  turning  point,  we  obtain 


sin  '  bx,  -f- 


b_ 

<«r+  byr- 


(2  Sin  [6a  +  0,(1  -  bJx/)]i/J)  <9) 


For  a  hard  sphere,  a2  is  large  compared  to  />  over 
most  of  the  scattering  range,  except  near  forward 
scattering  where  b  achieves  its  maximum  value  of 
I  /x2.  We  may  then  neglect  the  second  arc  sine  term 
in  (9)  and  write 


X1  6  .  1  ,  .  w  b 

—  =  sin  bx,  +  ,  (10) 

2  2  a. 

Solving  for  b  and  taking  b/a-/'-  as  small  compared 
to  unity,  we  obtain 


.  _  _  cos  (0/2) _ 

*„[!  +  (ir/2MW),/lf  sin  (0'2)| 

which  yields  the  cross  section 


(ID 


*  I  +  (x/2)  csec  (6/2)(a,x,!)  1 
x,u  [I  +  (r/2)  sin  (6/2)(aW)  V'J)'  ' 


An  analysis  of  (9)  and  (11)  shows  that  the  cross- 
section  formula  given  by  (12)  is  valid  for  hard 


!S*> 

spheres  from  the  backscattcr  region  to  about  30°  in 
the  forward-scatter  region.  Therefore,  to  complete 
our  calculation,  we  require  the  cross  section  for 
small  scattering  angles  in  the  forward  direction. 

Because  we  arc  interested  in  small  values  of  0  and 
hence  values  of  .r  where  the  electron  density  is  small, 
we  approximate  /e  by  the  linear  function 

a’  —  I  -  u,i(.v  —  .*„)  (13) 

where  .v„  is  a  cutoff  position  at  which  the  electron 
density  is  negligibly  small.  The  coefficient  an  is  de¬ 
termined  from  the  slope  of  the  a2  versus  x  curve  at 
v„;  that  is, 


o„  —(2,  x„)  exp  (,tor  2  -  ,r0  a)  (14) 
The  scattering  integral  is  then 


v  —  0  .  , 

—  =  sin  bx„ 

+  bLu  (,5) 

For  large  impact  parameter,  xt  is  close  to  x0\  we 
may  therefore  approximate  1  -  b* x*  by  2(1  -  fur) 
to  give 


sin  1  bx„ 


J„  ((2  +  O  -  (2b  +  a0)x}W! 
=  sin  1  bxn  +  [26  (2*  +  «„)][2(1  -  A.r„)],/2 
The  term  bxn  is  close  to  unity;  therefore,  we  can  write 


sin  1  bx0  S  *  -  [2(1  -  for.)]"* 

(16) 

$  =  (l  +  — )  [2(1  -  fur„)l''-‘ 

2  \  OnxJ 

which  yields  the  cross  section 


It  is  interesting  to  note  that  no  change  is  produced  in 
(16)  if  we  use  in  addition  to  the  linear  term  in  ( 13 ) 
a  quadratic  term  bu(x  -  .r0)2.  This  result  is  not  un¬ 
expected;  the  integrand  is  linear  in  the  small  quantity 
I  -  bx  and,  therefore,  terms  of  higher  order  than  the 
first  in  the  small  quantity  (x  -  x0)  arc  negligible. 

Because  the  cross  section  becomes  infinite  for 
zero-scattering  angle,  it  is  convenient  to  adopt  a  cut¬ 
off  value  for  the  impact  parameter  so  that  a  large  but 
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finite  crou  section  is  obtained.  A  reasonable  require* 
ment  is  to  take  the  electron  density  N0/Net  at  x0  be¬ 
tween  l(h*  and  10"*.  Using  (5)  and  (14)  we  express 
the  cross  section  at  0  =  0  by 

*">  -  ('  +  »5k)  (l8> 

where  x0  is  determined  by 

*<r2  =  Arcr“J  +  log  (A.r/A'o)  (19) 

3.2.  Intermediate  case.  If  the  sphere  is  mod¬ 
erately  hard,  the  critical  value  xct  will  be  somewhat 
greater  than  the  point  of  inflection  which  will  affect 
the  right  side  of  the  curve.  We  take  as  typical  ex¬ 
amples  r„  =  1  and  rct  =  1/2  for  which  plots  of  /** 
against  x  are  shown  in  Figures  3  and  4.  For  the  case 
r„  —  1  (Figure  3)  we  see  that,  aside  from  the  bend¬ 
ing  in  the  region  n  s*  1,  the  curve  is  almost  a  straight 
line  over  the  remaining  region.  Although  it  is  possible 
to  use  a  parabola  over  the  initial  portion  of  the  curve 
and  a  straight  line  over  the  remainder,  the  resulting 
formulas  become  sufficiently  complicated  to  make 
analysis  and  numerical  calculation  difficult.  We  shall, 
therefore,  use  a  straight  line  starting  at  x  -  0.54, 
M  =  1  to  jr,.r  =  1.0,  fKcr  =  0.  The  initial  portion  of 
the  curve  from  x  =  0.54  to  x  =  0.6,  while  not  ac¬ 
curate  should  not  affect  the  remainder  very  much 
since  it  is  a  small  part  of  the  integral  in  (4)  over 
most  of  the  range.  Furthermore,  because  it  is  above 
the  exact  curve  in  the  initial  region,  this  should  com¬ 
pensate  somewhat  for  the  straight  line  being,  on  the 


Fig.  3.  Plot  of  m*  for  moderately  hard  sphere  (r„  =  t ) 
and  approximation  by  single  straight  line. 


Fig.  4.  Plot  of  ti2  for  moderately  soft  sphere  (r„  =  }) 
and  approximation  by  two  straight  lines. 


average,  slightly  below  the  exact  curve  in  the  re¬ 
maining  region.  Our  straight  line  can  be  written  as 


a  —  l  rSx, 

ti1  =  1  —  «, (at  —  .*i)  *  2:  At, 

for  which  the  scattering  integral  yields 


(20) 


sin  1  Ajc,  +  ^ 


—  sin 


fl|  -j-  2ft2Af| 


la/  +  46  (1  +  a,*,)]1 


(21) 


The  arc  sine  terms  may  be  combined  in  a  simple 
manner  and  the  parameter  b  solved  for  explicitly  as 
a  function  of  scattering  angle.  We  obtain 


b  =  x,'1  cos  (0/2) 


Tl  +  —  (l  +— )  sin2 1/2 
L  a,*,  \  a,x,/  2J 

which  yields  the  cross  section 

1  +  2/ Q|AC| _ 


(22) 


v(0) 


JL 

~~  2 

*\  L 


I  +  4/fl,*,(l  +  l/a.-Osin2  (0/2) 


I 


(23) 


An  analysis  of  (19)  shows  that  when  the  turning 
point  =  0.6,  the  corresponding  angle  is  less  than 
30°;  thus  the  cross-section  formula  can  be  used  to  at 
least  30°.  The  cross  section  near  zero  scattering  can 
then  be  obtained  from  (17). 

As  the  critical  radius  decreases  to  0.5,  Figure  4 
shows  that  the  right-hand  portion  of  the  ^  -  x  curve 
becomes  concave  upward  so  that  the  curve  may  no 
longer  be  fitted  accurately  by  a  single  straight  line.  It 
is  possible,  however,  to  obtain  good  accuracy  by  the 
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use  of  two  straight  lines  of  different  slopes  meeting  at 
t  =  1.4.  The  cross  section  corresponding  to  turning 
points  within  the  region  of  the  first  straight  line  may 
be  obtained  from  (23).  When  the  turning  point  lies 
beyond  the  region  of  the  first  straight  line,  the  contri¬ 
bution  to  the  scattering  integral  from  the  second  line 
must  be  considered.  Because  of  the  discontinuity  in 
the  slope,  the  integral  loses  its  simple  structure  and 
can  no  longer  be  solved  in  simple  analytic  form  for 
arbitrary  b.  However,  for  small  b  (large  scattering 
angles),  the  integral  assumes  a  simple  form  and  can 
be  conveniently  solved.  For  turning  points  beyond  x2 
the  scattering  integral  is 


,  a.  +  2  b\ _ 

bxx  +  sin  [a,  +  4A„(|  +  ij  x  jj-i 


_ a,  +  26*, 

la,1  +  4b\\  +  a,*,)]"* 


,  r  .  _ a,  +  2b' Xj 

+  2  ,,n  [a1  +  46  W  +  atXi)\'u 


(24) 


For  small  b,  the  arguments  of  the  arc  sine  terms  may 
be  expanded  to  first  order  in  b 2  to  yield 


it _ b  ■  -i  ,  ,  ~i( ,  2 b 

— —  =  Mn  bx,  +  sin  y  ~  V?  J 

~  Sin"‘  (*  ~  +  I  “  Sin"  (‘  -  af3  ^0  <25> 


where  we  have  used 


»i‘  -  1  -  a,(jra  -  jc,) 

Expanding  the  arc  sine  terms  for  small  b,  we  obtain 


*  -  6  r 

-‘V' 


+  —  (I  — 

a,x,  a a 


which  yields  the  coss  section 


.r,  sin  6 


*> 

1  +—  (I  - 
L  a,.v, 


valid  for  large  scattering  angles. 

3.3.  Soft  sphere.  When  the  critical  radius  is 
small,  the  initial  moderate  decrease  in  ^  is  followed 
by  a  very  slow  decrease,  as  shown  in  Figure  5  for 
rcr  =  1/8.  It  becomes  difficult  and  inaccurate,  there¬ 
fore,  to  use  a  combination  of  two  simple  curves. 
However,  we  have  found  it  feasible  to  use  a  straight 
line  for  the  initial  portion  of  the  curve,  and  an 
approximate  method  of  integration  for  the  latter 
portion. 

A  reasonable  fit  up  to  x  =  1.4  may  be  obtained 
with  a  straight  line  from  *1  =  0.5,  /*2  =  1.0,  to  x2  — 
1.4.  The  cross  section  may  then  be  obtained  from 
(22)  and  (23)  as  long  as  the  parameter  b  does  not 
go  below  the  value  corresponding  to  x2  as  the  turn¬ 
ing  point,  that  is,  b  =  ii2/x2.  The  calculation  then 
yields  a  maximum  scattering  angle  of  about  69°. 

For  values  of  x  beyond  x2  we  split  the  integral  at 
a  suitable  point  .r:t  where  xs  is  not  far  from  x2,  but 


o  I  *  j  4  s  •  r  • 


X 

Fig.  5.  Plot  of  m2  for  very  soft  sphere  (r,,  =  V» )  and  approximation  by  straight 
line  for  initial  position  of  curve. 
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such  that  .r3J  is  somewhat  larger  than  unity  and  write  which  yields  the  cross  section 


/  =  /,  +  /..  (28a) 

where 

« 

'•  -  £ ;  ;■---??)»*  ,Me> 

The  integral  /,  may  be  evaluated  approximately  by 
noting  that  n  has  not  changed  very  much  in  the  in¬ 
terval  and,  therefore,  may  be  taken  as  constant.  We 
obtain 

.  .  i  f».v»  .  i  bx, 

I,  =  sin  T-r  —  sin  (29) 

(p>  <m) 

where  (ji)  is  an  average  value  in  the  interval.  To  in¬ 
tegrate  h .  we  assume  .tv  is  sufficiently  greater  than 
unity  so  that  a  scries  expansion  of  ^  may  be  used. 
We  then  have 

V  =  1  -  exp  (,x,,  2  -  .x'2)  (30a) 

£  -t'2  -  -v„"s  (30b) 

The  integral  /2  now  takes  the  form 

,  .  f" _ (3n 

2  J„  (x..,‘  -  X2  -tb'X'SxY*  ’ 

and  may  be  integrated  to  yield 

,  _  *  1  .  -i  1  +  2 b'x„2x, 

/2~4~2S,n  (1  +  4.xrrV)‘72  (  32) 

For  values  of  9  greater  than  and  near  ir/2,  the  b: 
term  in  the  numerator  is  comparable  to  unity  whereas 
the  b 5  term  in  the  denominator  is  much  greater  than 
unity.  It  is  not  until  6  is  in  the  range  near  ir  that  the 
latter  term  becomes  comparable  to  unity.  We  find  it 
convenient,  therefore,  to  calculate  near  v/2  and  near 
w.  For  6  near  v/2,  the  arc  sine  terms  in  /t  and  l3  are 
small  compared  to  unity,  and  we  may  write  our  scat¬ 
tering  integral  as 

t/2  -  6  =  (A  -  a-,2)  b  -  (2j (33) 

where  A  =  2[xt  +  (jr3  -  *2)/(/i)].  Solving  for  b,  we 
obtain 

[M’ 

+  4<<l^l>]'-}  (J4) 


sin  ex„\A  -  .X:,2)2 

{(!-»)  +  [M’  +  ’-f/T} 

[('-*)'  + 2  ~2‘ ‘-T  051 

In  terms  of  b,  the  cross  section  has  the  simpler  form 

'<*>  -  ;?ne  I  +  2  ,36) 

For  large  scattering  angles,  the  argument  of  the 
arc  sine  term  in  the  I3  integral  is  no  longer  small 
compared  to  unity;  thus  the  preceding  simplification 
cannot  be  used.  We  therefore  write  the  scattering  in¬ 
tegral  in  the  form 

«  4L  j  ,1  +  2 b*xr,W 

■  -  »  -  Ab  +  COS  ^-.^,7,  im 
and  take  the  cosine  of  both  sides  to  yield 

”u,,o  .KV,“ 

-a,i,l*'(}TSi) 

Since  Ab  is  a  small  quantity,  we  may  expand  cos 
(Ab)  and  b  sin  (Ab)  to  first  order  in  b-  and  solve 
for  b  to  obtain 

b  =  sin  (it  —  9)/2xf/[B  +  cos2  Or  —  fl)]l/2  (38) 

where 

n  -  ~A,  (t  -  -sV  -  , 

■X,  ,  \  .xtr  /  ,v,t  4a,  , 

The  cross  section  is,  accordingly, 

it  (I  +  D)  cos  (t  -  9) 

o(9)  =  — 4-  (39) 

.xrr  [/?  +  cos  (»  -  0)] 

For  backscattcr,  we  obtain  the  interesting  limiting 
form  (neglecting  the  small  quantity  B) 

o(8)  =  ir  .x,/  (40a) 

r,.n  =  r,,1  (40b) 

where  r,.tt  is  the  radius  of  a  conducting  sphere  having 
the  same  cross  section  as  the  Gaussian  sphere. 

3.4.  Underdense  sphere.  When  the  maximum 
electron  density  of  the  plasma  is  below  the  critical 
frequency,  the  index  of  refraction  has  a  minimum 
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value  greater  than  zero.  The  general  shape  of  the 
curve  is  similar  to  a  soft  sphere,  with  the  minimum 
value  approached  asymptotically  for  large  .r.  In  view 
of  the  tediousness  of  our  method  for  very  soft  spheres, 
we  use  a  perturbation  technique  in  which  we  assume 
in  first  approximation  that  the  deflection  of  the  ray 
path  may  be  neglected.  For  a  spherical  plasma,  the 
relation  between  index  of  refraction  and  the  angle  a 
between  the  radius  and  the  ray  direction  is  given  by 
Bouger’s  rule  [KV/jo,  1964] 

pr  sin  a  =  b  (41) 

The  variation  in  «  due  to  change  in  n  is 


a  relation  similar  to  that  obtained  in  classical  mechan¬ 
ics  [ Bohm .  1951]  if  we  take  p  as  close  to  unity.  For 
our  Gaussian  sphere  p  I  -  pe~’’  where  p  - 
AL,»  'N„.  yielding 

<«> 

J  »  I  —  pe  e 

where  u  -  z/h.  Since  p  is  less  than  unity,  we  may 
expand  the  denominator  to  effect  the  integration  and 
obtain 

6  =  rutcb[\  +  -*m  +  -4  +  •■■)  (45) 


da 


tan  0  d(p') 
— — T~  dr 

2p  clr 


(42) 


where  we  omit  the  change  in  «  due  to  the  explicit  r 
term  since  this  gives  its  variation  due  to  the  angular 
change  in  r  (note  that  for  fixed  p  and  hence  no  de¬ 
flection,  the  variation  in  r  gives  precisely  the  change 
in  a).  Since  da  is  referred  to  a  fixed  position  of  r,  it 
gives  the  deflection  of  the  ray  at  this  point.  The  total 
deflection  of  the  ray  is  then  obtained  by  integration 
over  the  ray  path.  If  we  assume  the  path  is  unper¬ 
turbed  then  tan  a  =  b/z,  r-  -  b-  +  zs,  and  the  in¬ 
tegral  can  be  cast  in  the  form 

e  =  (43) 
2  J.m  n  dr  r 


where  c  pe'".  which  yields  the  cross  section 

a(0)  =  A^-P  <L[~(|  -  2bJ)  +  (1  -  4ft2)  TT/j 
sin  0  l  L  4 

+  (1  —  6ft’)  jpii  +  ]}  (^) 

For  underdensc  spheres,  the  scattering  angle  vanishes 
at  b  =  0  as  well  as  at  large  b  so  that  b  is  a  double 
valued  function  of  6\  accordingly,  the  contributions 
from  the  two  values  of  b  must  be  added  numerically 
to  obtain  the  total  cross  section.  The  scattering  in 
the  forward  direction  may  be  calculated  by  the  pro¬ 
cedure  for  overdense  spheres.  However,  since  P  is 
less  than  unity,  we  express  (*9)  directly  in  terms  of 
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Fig.  6.  Variation  of  impact  parameter  with  scattering  angle  for  overdense  spheres 
for  levertl  values  of  critical  radius. 
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Fig.  7.  Variation  of  effective  radius  with  scattering  angle 
for  several  values  of  critical  radius. 


p  by 

jro'3  =  \n(Sc,/No)  ~  In  P"  (47) 

4.  RESULTS  AND  DISCUSSION 

The  results  obtained  by  the  numerical  procedure 
for  overdense  spheres  are  presented  in  Figures  6  and 
7,  where  we  have  plotted  the  variation  of  impact 


Fig.  8.  Plot  of  ratio  of  effective  radius  to  critical  radius 
against  scattering  angle  for  overdense  spheres  with  critical 
radius  as  the  parameter. 


parameter  and  cross  section  with  scattering  angle, 
with  the  cross  section  expressed  in  terms  of  rpff.  The 
results  obtained  by  the  analytic  procedure  were  gen¬ 
erally  very  close  to  those  obtained  numerically.  To 
indicate  the  agreement  obtained,  we  have  also  plotted 
in  Figure  7  the  analytic  solutions  for  re„  for  a  hard 


Fig.  9.  Variation  of  impact  parameter  with  scattering  an¬ 
gle  for  nnderdense  spheres  for  several  values  of  density  ratio. 
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sphere  (r„  =  2)  and  for  a  very  soft  sphere  (r„  = 

1/8). 

Figure  6  indicates  that  for  hard  spheres,  except 
for  the  forward-scattering  region,  the  impact  param¬ 
eter  is  almost  a  linear  function  of  scattering  angle. 
As  the  critical  radius  becomes  small,  the  impact  pa¬ 
rameter  decreases  rapidly  in  the  forward-scattering 
region,  resulting  in  very  small  values  for  all  scatter¬ 
ing  angles  greater  than  90°. 

Figure  7  shows  that  hard  spheres  have  an  almost 
constant  cross  section  over  a  fairly  wide  range  of 
backscatfering  angles,  and  rt.„  is  only  slightly  less 
than  r„.  For  r(t  =  0,  the  cross  section  becomes  van¬ 
ishingly  small  from  0  =  90°  to  180°.  We  have  in¬ 
dicated  this  result  by  a  vertical  line  at  90°.  The  small 
values  of  rM  for  a  soft  sphere  in  the  backscatter  re¬ 
gion  are  due  to  the  great  depth  of  penetration,  and 
consequently  small  value  of  the  impact  parameter  b 
required  for  a  given  scattering  angle.  The  rate  of 
change  of  impact  parameter  with  scattering  angle  is 
accordingly  very  small  in  the  backscatter  region,  re¬ 
sulting  in  low  cross  sections.  We  thus  have  a  transi¬ 
tion  region  in  which  reM  decreases  very  rapidly  until 
the  scattering  angle  becomes  quite  large.  This  result 
may  also  be  seen  from  the  form  of  the  /r  versus  x 
curve  for  a  very  soft  sphere  (Figure  5)  and  the  scat¬ 
tering  integral.  As  the  turning  point  gets  closer  to  the 
center  of  the  sphere,  x,  becomes  larger  and  the  inte¬ 
gral  starts  to  increase  rapidly  because  of  the  large 
contribution  from  the  region  of  small  n  values;  the 
valu  j  of  b  then  begins  to  decrease  rapidly  with  a 
sharp  drop  in  cross  section.  For  large  values  of  the 
turning  point,  however,  the  integral  is  quite  large  and 
further  increases  in  the  scattering  angle  have  very 
little  effect  upon  db/d6,  resulting  in  a  slowly  decreas¬ 
ing  cross  section. 

In  the  region  of  forward  scatter,  we  note  that  the 
curves  for  different  r„  coalesce,  tlr  process  occur¬ 
ring  at  smaller  angles  for  larger  r„.  This  result  occurs 
because  the  forward-scattered  rays  penetrate  the 
outer  edge  of  the  cloud  where  the  electron  density 
and  its  gradient  are  small,  and  do  not  differ  signifi¬ 
cantly  for  different  critical  radii.  The  scattering  there¬ 
fore  depends  principally  upon  the  impact  parameter 
and  very  little  upon  the  critical  radius.  For  large  val¬ 
ues  of  r<r,  however,  the  electron  density  is  affected 
to  some  extent  by  the  critical  radius  unless  the  im¬ 
pact  parameter  is  taken  sufficiently  large.  (See  equa¬ 
tion  19,  where  we  note  that  large  r„  may  have  some 
effect  unless  the  cutoff  density  N0/Nt ,  is  taken  suffi¬ 
ciently  small.)  The  curves  for  large  rt„  therefore, 


coalesce  more  slowly  than  do  the  smaller  r„  curves. 

Since,  in  backscatter,  the  effective  radius  reff  is 
very  close  to  r„  for  large  rcr  and  decreases  as  rer: 
for  small  rcr  (see  equation  40),  it  is  of  interest  to  see 
whether  a  simple  exponential  law  can  be  used  to  rep¬ 
resent  the  backscatter  results.  A  plot  of  log  (1  — 
»VffAer)  against  r„  shows  that  the  backscatter  results 
may  be  fitted  with  good  accuracy  by  a  straight  line 
whose  equation  is 

1  ~  r.,,/r„  =  e‘,'0,,,,  (48) 

This  reli  tion  gives  the  proper  dependence  of  r,ff  upon 
r„  for  small  rf,  and  satisfies  the  requirement  that 
rttt/r„  is  close  to  unity  for  large  r„.  To  show  the  re¬ 
lations  between  rttt  and  r„  for  all  scattering  angles, 
we  have  plotted  in  Figure  8  the  ratio  rttt/rCT  as  a 
function  of  scattering  angle  for  several  values  of  rct. 


Fig.  10.  Variation  of  effective  radius  with  scattering  angle 
for  underdense  spheres  for  several  values  of  density  ratio. 
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The  curves  cross  near  8  -  105  '  and  r,.„  r,.,  indi¬ 

cating  that,  in  this  scattering  region,  all  Gaussian 
spheres  may  be  replaced  by  a  conducting  sphere  of 
radius  equal  to  the  critical  radius. 

The  variation  of  impact  parameter  with  scattering 
angle  for  underdense  spheres  is  shown  in  Figure  9 
where,  in  contradistinction  to  the  overdense  case, 
small  deflections  arc  obtained  at  small  as  well  as  at 
large  values  of  the  impact  parameter.  The  results  ob¬ 
tained  by  the  perturbation  procedure  are  in  fair  agree¬ 
ment  with  the  more  exact  calculations  at  small  values 
of  the  density  ratio  P,  the  agreement  lessening  as  p 
increases.  The  more  exact  calculations  show  a  steady 
decrease  of  b  at  maximum  deflection  with  increasing 
p,  whereas  the  corresponding  value  of  '<>  by  the  ana¬ 
lytic  method  is  relatively  insensitive  to  an  increase  in 
p.  The  two  sets  of  curves  are  therefore  considerably 
displaced  at  higher  values  of  p,  with  consequently 
inaccurate  cross  sections.  Attempts  to  improve  the 
perturbation  results  by  an  iterative  procedure  in 
which  the  orbit  is  not  neglected  have  not  proved 
very  successful.  We  have  therefore  not  presented  the 
analytic  results  for  the  cross  sections. 

The  cross-section  results  are  presented  in  Figure 
10,  where  we  have  added  the  contributions  from  the 
two  values  of  b.  The  contribution  from  the  large 
value  of  b  is  generally  much  higher,  except  near  0,1UU. 
where  both  b  values  give  very  large  contributions. 
The  vertical  tines  in  Figure  10  indicate  the  points  of 


maximum  deflection  at  which  the  cross  section  is  in¬ 
finite.  As  anticipated,  the  curves  for  different  p  coa¬ 
lesce  in  the  forward-scattering  region,  with  merging 
occurring  at  lower  values  of  8  for  smullcr  p.  This 
may  he  easily  seen  by  noting  that  the  value  of  r„  de¬ 
pends  more  strongly  on  p  when  it  is  small  (sec  equa¬ 
tion  47 1. 

For  the  purpose  of  comparison  and  discussion  with 
respect  to  the  critical  density  region,  some  of  the  im¬ 
pact  parameter  results  tor  overdense  uud  underdense 
spheres  arc  presented  together  in  Figure  1 1 .  We  note 
that,  as  the  density  approaches  unity,  the  upper 
branch  for  the  underdense  region  merges  with  the 
forward-scatter  portion  of  the  overdense  curve.  In 
addition,  the  lower  branch  of  the  underdense  curse 
goes  to  h  =  0  with  a  maximum  deflection  of  90'  . 
while  the  backscatter  portion  of  the  overdense  curve 
also  goes  to  b  --  0  for  the  region  90°  to  180°. 

We  thus  have  a  discontinuity  in  maximum  scatter¬ 
ing  angle  at  the  critical  density,  with  90°  obtained 
when  the  approach  is  from  underdense  values  and 
1 80'  when  it  is  from  the  overdense  side. 

The  cross  section  approaches  zero  in  the  back- 
scatter  region  as  the  density  becomes  critical  from 
the  overdense  side  (see  Figure  7)  whereas  because 
of  the  infinite  slope  in  the  impact  parameter  curve, 
it  becomes  infinite  at  90°  from  the  underdense  side 
(Figure  9). 

These  peculiarities  in  the  scattering  near  critical 


Fig.  1 1 .  Plot  o(  impact  parameter  against  scattering  angle  for  underdense  and  over- 
dense  spheres  showing  limiting  forms  at  critical  density. 
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density  are  due  to  the  use  of  a  geometric  optics  model 
for  the  calculations.  The  nature  of  the  scattering  in 
this  region  may  be  exhibited  in  terms  of  a  simple 
model  for  the  electron  density  which  has  the  proper 
asymptotic  form  near  the  center  of  the  cloud  ( large 
,r ) .  We  write 


n‘  =  1  —  pa‘x‘  0  S  j  S  x, 

i 

m'  =  I  —  p  4-  “r  -v  =  x. 


(49) 


where  a  simple  parabolic  fit  has  been  chosen  for 
small  .v  and  an  inverse  square  behavior  for  large  i. 
The  coefficients  a  and  c  may  be  evaluated  in  terms 
of  a  specific  density  distribution. 

The  scattering  integral  now  reads 


2b 

(up  +  by 


/  *  1 

(a  p  b  )  x% 


v  .  2 b'x,‘  —  (I  —  p) 

+  2*sm  [Ti  -  p^pc-vr 


(50) 


For  b  close  to  zero  the  first  term  is  negligible  indicat¬ 


ing,  as  expected,  that  the  scattering 
parameters  is  determined  by  the  inti 
center  of  the  sphere.  For  b  =  Owct 

7T  —  9  =  (r/2)  —  sin  '  [(I  —  p)/| 

which  yields  9  —  0°  for  underdense 
180°  for  the  overdensc  case.  At  tlu 
the  limiting  value  as  b  approaches 
addition,  the  model  gives  the  requiret 
of  b  versus  0  near  9  =  0  ( underdens 
18(T  lovcrdense). 
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